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The closed form solution for the geodesies of classical particles in SdS space 
is obtained in terms of hyperelliptic modular functions and multiple hypergeo- 
metric functions. The closed form solution for the five roots of the fifth degree 
polynomial is found giving the branch places on the genus two Riemann surface. 
'The Inversion Problem', for the genus two hyperelliptic integral, is solved in 
a closed form. The solution is shown to reduce to elliptic functions when the 
cosmological constant is zero. Current observational data is in favor of a cosmo- 
logical constant model 1 . This solution is important in astrophysical applications 
of measuring the cosmological constant. 

I. Introduction 

To calculate the geodesies in the centrally symmetric gravitational field of a 
Schwarzschild black hole in a universe with a cosmological constant(A 7^ 0), 
i.e. Schwarzschild de-Sitter space(SdS space), the Hamilton- Jacobi method is 
applied to the metric 2-5 

ds 2 = g tt dt 2 + g rr dr 2 + g^dep 2 + geedO 2 (1) 
g tt = (l--erAr 2 — )c 2 ,g rr = -(1-icrAr 2 -^) -1 ,^ = -r 2 sin 2 9,g ge = -r 2 

O V o V 

(2) 

Where a = ±1 depending on whethor or not the universe is de-Sitter or Anti- 
de-Sitter, respectively. A is the cosmological constant and was first introduced 
by Einstein 6 and r g was introduced by Schwarzschild. The functional form of 
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the geodesies when r g — have been studied, they are given in terms of genus 
zero circular functions. When A = the geodesies are given by genus one 
elliptic functions 7 . In this paper , the functional form of the geodesies when A 
and r g are both non-zero are found, in terms of genus two hyperelliptic thcta 
functions 8 ~ 9 and multiple hypergeometric functions 10 . 

II. Derivation of the Hyperelliptic 0-Integral 

Let a particle of mass m be constrained to move in the equatorial plane, 8 = 0, 
and let the action be given by 11 , 

S = -Et + L<f> + S r (3) 

where E is the total energy of the particle and L is it's angular momentum. By 
the conservation of four momentum, 

• • dS dS 09 . . 

° V d*W =m<? (4) 

we can insert S from (3) and g 1 ^ from (2) into (4) and obtain, 

Sr ~L (i - \^ - 5) • ™ 



Differentiating S r w.r.t. L and using 

, dS r 



const = (6) 
oL 



-f 



Ldr 



(7) 



r 2 (f - (1 - iAar2 - ^)(^ + m 2 c 2 ))5 

This integral is a hyperelliptic integral of genus, g = 2. The problem is to invert 
this integral, i.e. express the upper limit r as a function of 0, 12-14 and as a 
function of r, this is known as the 'Inversion Problem'. The inversion problem 
was first conceived by Euler 15 . 

III. The Solution to the Inversion Problem for the 
Hyperelliptic Integral of Arbitrary Genus 

The following will be an account given by H.F. Baker 8-9 and H. Exton 10 on the 
application of theta functions and multiple hypergeometric functions to obtain 
the closed form solution to the arbitrary genus hyperelliptic integral, except for 
the branch places on the Riemann surface. This requires the solution to the 
n th d e g ree polynomial, for its n roots. The general solution to the inversion 
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problem is given by solving g of the 2g + 1 equations, for the g variable places 

X\ , . . . , Xg , 

^± F J ) = A(b)(b- Xl )...(b-x g ) (8) 
The notation Baker uses for the generalized theta function is, 

0(u|u M ) = ^H^fW) - 0(u; \m, im') - e a " 2 e( W ; ^m, ^m') (9) 

where a is an arbitrary g x g symmetrical matrix, since putting (9) into (8) 
gives, 

where m, m', k and fc' are integers. Riemann's theta functions are defined as, 
e( u; I m , im') = s e 2/ "' (Tl+ i m ' )+(,(n+ i m ' )2+i ' rm(n+ i m ' ) (11) 

with, X = Y^inv=-oo ■ ■ ■ En 9 ™-oo' h is & g x g matrix, in general asymmetrical. 
b is a # x g symmetrical matrix. In (10), u denotes the g quantities, 

u? 1 ' 01 + . . . + u^ a <> = Ui (12) 

with i = 1, . . . , g. This is known as Abel's Theorem. Riemann's Normal Integral 
of the first kind is defined as, 

1 -J a y { } 

where 

{x, l) i;S _i = A^g^x 9 - 1 + ... + A h0 (14) 

and the Abel coefficients A^ 9 _i, . . . , Ai,o are constants determined experimen- 
tally. The denominator in (13) is given by, 

y 2 = 4(x -ai)...(x- a g )(x -a)...(x- c g )(x - c) (15) 

where ai, . . . , a g , Ci, . . . , c are branch places and x is a variable on the one com- 
plex dimensional Riemann surface of genus, g. Now since the sum of terms is 
finite in (13), then by combining (13)(14)and(15) an arbitrary term in u^' a may 
be transformed to Extonian form 10 , using a Mobius transformation, 

j Z t^ a - l \l - _ ... (1 - tz 2g ) { -^dt 

Jo 

= ^-i^ 29+1) (a, -,1 + a-c; a + l;zi,...,z 2g ,z) (16) 

where Re(a) and Re(c — a) are greater than zero, a and c are appropriately 
chosen constants. The variables z\,...,z n depend on the roots of y and the 
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lower limit. While z depends on the upper limit of the hypcrclliptic integral. 
The condition that z\ , . . . , z 2g and z each have their absolute value less than one 
ensures that the Lauricella function F^ 9+1 ^ converges, making the variables in 
Abel's theorem well defined. The convergence takes place in a unit sphere in 
2g + 1 hyperspace. The Lauricella function has the expansion, 

F% 9+1) (a, i, . . . , 1 + a - c; a + 1; zi, . . . , z 2g , z) 

(o, mi + . . . + m 2g +i){\, mi) ... (1 + a - c, m 2g+1 )z] ni . . . z m ^ 
. f-^ (a + l,mi + . . . + m 2g+ i)mi! . . .m 2g+ i! 

(m)=— oo 

|zi|<l,...,|^ 2ff |<l,|z|<l (18) 
where (m) = mi . . . m 2g and have the form from the Mobius Transformation, 

Zi = (19) 

n-a 

where a is the lower limit of the hyperelliptic intcgral(13), and 



Now let, 



x — a 

z = (20) 

n-a 



^ m ,m> = miu) r ,i + • • • + m g u) r ,g + m[uj' r l + . . . + m' g J rg (21) 

where the uu have the representation, 

w = C x F^-Via, i, . . . , i; c; zi, . . . , z 2g _i) (22) 

and where m and ml are equal to +1 when integrating across a period loop 
right to left, and —1 when crossing from left to right and when not crossing 
a period loop. C is a constant that depends on the branch places and Abel's 
coefficients, lu is convergent in a unit sphere in 2g ~ 1 hyperspace. a and c are 
the appropriate constants. Now define an integral of the first kind via, 

Triv^ = //,.,„;-• + . . . + h ri9 u x g ' a (23) 

with r = 1, . . . , g and where, 

2hu = m (24) 
2hu' = irir (25) 
where Riemann's condition for hyperelliptic integrals is imposed 16-18 . 

T = T* (26) 
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With the above equations, Riemann's theta function becomes, 

Q(v; —m, -to') = ^ e 2viv ( n +i m ')+ i7TT ( n +i m ') 2 + m ( n +i m ') (27) 

where the quadratic forms are given by, 

v(n + ^m') = vn + ^m'v = vini + . . . + v g n g + \^\m! x + . . . + ^v g m' g (28) 

and similarly, 

r(n + -to') = rn 2 + mm 1 + -rm' 2 
2 4 

= (Ti.inj + 2ri i2 nin2 + . . . + T g , g n g 2 ) 

a a 

+ X! T r,sn r m' s 



5=1 r— 1 

+ J(n,in»i 2 + 2ti, 2 to' 1 to 2 + . . . + t s , s to^ 2 ) (29) 

It has been proven by Riemann that for the integral of the first kind, when it 
has the given period scheme, then the imaginary part of, 

ri,im 2 + 2Ti i2 nin 2 + . . . + T g ^ g n 2 (30) 

is positive for all integer values of n\ and n 2 , with the exception of n\ — n 2 = 0. 
Therefore the modulus of e 47rr " is less than unity and the function(27) converges 
for all values of its argument, v. The last thing to define in (8), is the constant 
A(b), it is defined as, 



A (b) = (e^(x-a 1 )...(x-a g ){x-c 1 )...(x-c g )(x-c)\ x=b ) i (31) 



where e is +1 when u > a is an odd half period and —1 when it is an even half 
period. 



IV. The Complete Solution to the Inversion Prob- 
lem of the Genus Two Hyperelliptic Integral 

The solution to the inversion of the genus two hyperelliptic integral consists of 
two parts. The first part is to express the upper limits, x\ and x 2 , in terms of 
the anglers. u\ and it 2 ; m Abel's theorem, 

«i I1 ' ai +tti Ia,ea =ui (32) 

U2 I1,0l +U2 Xa, ° 2 =«2 (33) 
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and secondly, express the angles u\ and u 2 in terms of x\ and x 2 , so that u\, 
u 2 , x\ and x 2 become well defined coordinates on Rummer's quartic surface 19 . 
The solution to the first part is given by 20 ~ 22 ' 8 ~ 9 , with g = 2, 



6 2 (^i, v 2 ; \m 1 , \m 2 , \m\, \m' 2 ) 



A(b)(b - Xl )(b - x 2 ) (34) 



6 2 (wi,w 2 ; \k\, \k 2 , \k' x , \k' 2 ) 
where by (23) we have, 

• x.a 7 x,a 1 7 x,a /or\ 

irivi = mv x — h\,\u 1 + h\, 2 u 2 (35) 

■ x.a t x,a , 1 x,a fr,r\ 

mv 2 — mv 2 — h 2 ,\u x + li 2}2 u 2 (3o) 
where the hij are the elements of the following matrix, 

m I lo 2 . 2 -Wi >2 



h 

A V _w 2,l W14 

and 

A = det(w) ^ (37) 

Now the solution to the second part of the inversion problem is facilitated by 
first choosing, 

x 2 = a 2 (38) 

by (13), we have 

J ax y 

u 2 =u^= r (4»,° + ^,i*)<fe (40) 

J a\ V 



where, 



y 2 = 4x 5 + A 4 x 4 + A 3 x 3 + \ 2 x 2 + Aia; + A (41) 
which can be rewritten as, 

y 2 = 4(x - ai)(x - a 2 )(x - Ci)(x - c 2 )(x - c) (42) 

Now let (39) and (40) be written as, 

Uh „y = f X (A K0 + A K1 t) dt 

Jy v /4 ( i - r n ) (t - n 2 ) (t - r l3 ) (t - r u ) (t - r l5 ) 

where k = 1, 2. i\,i 2 , is, 14 and i$ £ H mod 5 where He [0, 1, . . . 00), so that, 

ro^i, r 2 ,r 3 ,r 4 <G (ai,a 2 ,c 2 ,ci,c) (44) 

The integral(43) may written as a product of factors, via the Mobius transfor- 
mation, 

z = (45) 
r»! -y 
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we have, 

u k x - v = B(A kfi +A kA y)x 

/ (1 - z)-*(l - zx i2 )~^(l - zx i3 )~i(l - zx i4 )~i(l - zx i5 )~^dz 
Jo 

+ B(r n - y)A kA x 

z(l - - zx t2 y^{\ - zxi 3 )~?(l - zx i4 )~^(l - zx i5 )~i(lz 



o 



(46) 



where 



B = B{y,r il ,r i2 ,ri 3 ,r u ,r i5 ) = — J- r- Ttl T (47) 

2i V (n 2 - y){r i3 - y)(r u - y){r l5 - y) 

and 

_ r n -y _ r h -y _ n.-y _ r il -y 
n 2 ~y n 3 -y r ti ~y r l5 - y 

The two integrals(46) arc of Extonian form, so that the intcgral(43)is given by, 

u k x ' y = phef(A kfi , A k ,i;x,y) + shef(A ki i;x,y) (49) 

where "phef ' stands for "primary hyper-elliptic function" and "shef ' stands for 
"secondary hyper-elliptic function". They are defined as, 

phef(A kfi ,A kA ;x,y) = B(A kfi +A k sy)zF D {5) (1, i, i, i, i, ^;2]x t27 Xi 3 ,x u ,x i5 , z) 

(50) 

shef(A kA ;x,y) = -^BA kA (r ix - y)z 2 F D (5) (2, 3; x l2 , x l3 , x u , x l5 , z) 

(51) 

Both (50) and (51) are defined via the twelve slotted, five variable, Lauricella 
D multiple hypergeometric function, Fd^ 9+1 \ with 2g + 1 = 5 Rcf. 10. These 
functions converge when the moduli of the five variable parameters Xi 2 , x is , Xi 4 , 
Xi 5 and z have their absolute values less than unity. The covergence takes place 
in a unit sphere in five dimensional hyperspace. In this case (50) and (51) may 
be expanded into their Lauricella series, i.e. 

phef(A kfi , A kA ;x, y) = B{A kfl + A kA y) x 

(l,mi + m 2 + m 3 + m 4 + m 5 )(±,m 1 )(±,m 2 )(±,m 3 )(±,m 4 ){^,rn 5 ) 
^ (2, mi + m 2 + TO3 + 77i4 + m^m\\m2\mz\m 4 \m^\ 

(m)=0 

v ic - m 2 T . m 3^,. rri4 r m 5 + l /M\ 

/\ ^0^2 ^-23 ~ V / 

shef(A kA ;x,y) = -BA kyl (r n - y) x 

(2, mi +m 2 +TO3 + m 4 + m 5 )(±, mi)(i, m 2 ){j, m 3 )(|, m 4 )(±, m 5 ) 
(3, mi + «2 + TO3 + m 4 + m 5 )mi!m 2 !m3!m4!m5! 

(m)=0 

v „,. ra x m 2 m s m 4 m 5 +2 /ro\ 

/\ <Ol 2 J-2 3 ^ V / 
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with 

\x l2 \ < 1; |ar i3 | < 1; |x u | < 1; \x ib \ < 1; \z\ < 1 (54) 
The quadrupoly periodic double theta function is given from (27) to (29), 

1111 oo oo 

e(wi,« 2 ;imi,^m 2 ,^m / i,^m / 2)= ^ ^ e ^ 1 (2„ 1 + m ' 1 )+^(2« 2 +m' 2 )) 

71 1 — — OO 7^2 — — OO 

x ^(27i 2 +m' 2 ) 2 ^/( 2n i+ m/ i) 2 r (2ni+m / i)(2n2+m / 2) j_|j(m 1 ni+m2n2)^(mim'i-m2m' 2 ) 

(55) 

where 

q = gi^i.i j g ' = e i^2,2 ) r = e |«n,2 (55) 
mi = p = p, ni2 = f> = A, rn'i = a = v , m' 2 = a' = p (57) 
so that a r in Forsyth's paper is equal to 

a r = (j)(°"'p'+pc r ) r (2n+CT')(2m+cr) ^gg-j 

where n = n\ and m = n 2 are the summation indices. Forsyth gives a review of 
Roscnhain's theory 22-23 of the fifteen ratios of the quadrupoly periodic theta 
functions. From the list given, the first ratio is selected and is equal to, 

9^,^1,1,0,1) / 0l -q 2 

e^v u v 2 ; 1,1,0,0) y (a, - Cl )( ai - c 2 )( ai - c) [ 1 j l ° yj 

where (31)(38)(55)and (59) have been used. The double theta function used in 
(59) is, 

+ CO + OC 

ni(v 1 (2n 1 )+v 2 (2n 2 + l)) 



*»= E E 

m = — oo n 2 — — oo 

x ? (2™ 2 +l) 2 9 /(2ni) 2 r (2n 1 )(2™ 2 + l) ( -_ 1 ^(n 1 + n2 ) i (-l) (gQ) 



where Forsyth's notation is adopted, 

fl 13 = e( t , 1 , V2 ;I,I,0,l) (61) 

and in the denominator of (59), is 

+ OC +OC 

# 12 = ^ ^2 e 7ri(t,l2 " 1+t ' 2(2 ™ 2)) (7 (2n2)2 (7' (2 " l)2 r (2 " l)(2 " 2) (-l)™ 1+ri2 (62) 

ni — — oo n 2 — — oo 

with, 

0i2 = e(«i,«2;^,O,O). (63) 

This completes the inversion of the genus two hypcrclliptic integral. The analytic 
continuation formulae are given in Exton. If one wants the solution with all four 
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coordinates x\, x 2 , Ui and u 2 in (32) and (33) one does not choose x 2 = a 2 - 
Proceed as above adding on the appropriate multiple hypergeometric functions 
and selecting another theta function. This was not done here to keep the solution 
as concise as possible while keeping its closed form (although all four coordinates 
may play a role in the geometric interpretation of the solution when applied to 
the problem at hand). Now in order to have a closed form solution we must find 
the roots of the quintic, n, r 2 , r 3 , r 4 and r 5 in a closed form. 

V. The Complete Solution to the Quintic 

A. The Tschirnhausian Transformation to Bring's Normal 
Form 

The initial Tschirnhausian transformation used is a generalization of Bring's 24-25 , 
but a simplification of Cayley's 26 , with a quartic substitution, 

Eql = x 4 + dx 3 + cx 2 + bx + a + y (64) 

to the general quintic equation, 

Eq2 = x 5 + tox 4 + nx 3 + px 2 + qx + r. (65) 

Then by the process of elimination between (64) and (65), the following 25 
equations are obtained, 

Mi 5 = 1, M M = d, Mi 3 = c, M12 = b, M u = a + y, M 25 = m-d 

M 2 4 = n — c, M 2 3 =P—b, M 2 2 = —y + q — a, M 2i = r, M 35 = n + dm-m 2 - c 
-^34 = p—b—mn+dn, M 33 = q—a—mp—y+dp, M 32 = r—mq+dq, M 3 i = dr—mr 
M45 = — cm— m 3 +b— p+dm 2 +2mn— dn, M 44 = a+dmn+y—dp—m 2 n—q+n 2 +mp—cn 
M42, = np — r+dmp — dq+mp — m 2 p—cp,M42 = np — r+dr — m 2 q — cq+dmq+nq, 

M41 = — m 2 r — cr + nr + dmr 

M55 = bm — 2mp + q — y + cn — 2dmn — a + 3m 2 n — cm 2 

+ dm 3 — n 2 — m 4 + dp 
M54 = —dmp + dm 2 n + cp + dq + 2mn 2 — cmn — m 3 n — 2np 

— mq + m 2 p + m 2 q — m 3 p + dr 

M53 = cq + 2mnp — dnp — p 2 + bp — nq + dm 2 p — mr — dmq 

— crap + m 2 q — m 3 p + r 

M 5 2 = bq — cmq — nr — dmr — dnq + cr + m 2 r — m 3 q + 2mnq 

— pq + dm 2 q — m 3 p + dr 

M 5 i = br — m 3 r — dnr + dm 2 r + 2mnr — pr — cmr 

(66) 
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These equations are then substituted into the equation, 



Mn 


M12 


M13 


M14 


Mib 


M 2 i 


M 22 


M 23 


M 24 


M 25 


M31 


M 32 


M 33 


M 34 


M 35 


M41 


M 42 


M 43 


M 44 


M 45 


M 51 


M 52 


M 53 


M 54 


M 55 



det(M) 



which generates a polynomial which will be reduced to Bring's equation(80). 
Let this equation be, 



V + P4V + PsV + PiV +Ay + B = 0. 



(67) 



Now take coefficients of the fourth, third, and second order terms in (67), re- 
spectively, and set them to zero. Solving, 



P4 = 



(68) 



for a, one has 



a = —dm 
5 



—bra dmn n 2 

5 5 5 



-q - -cm 
1 



2 4 4 2 3,,. 

H — cn mp H — m n H — dp — -m 

5 5 5 5 5 

Substituting a and Bring's substitution for b and c, i.e. 

b = ad + £ 
c = d + 1] 



(69) 

(70) 
(71) 

(72) 

with the third and second order coefficients in (67), p 3 and p 2 set to zero. P3 is 
a quadratic in d, so that in may be written as, 



into(67), which becomes 



y + Ay + B = 



P3 = P32d 2 + pud + p 30 = 0. 
Setting the coefficients, P3 2 , pzi and p 30 each to zero, gives 



(73) 



2 17 

P32 = (--m 2 + n)a 2 + (—m 2 n ■ 
5 5 



17 4 
—mp + -m 3 

5 5 

19 



2n z 
12 



13 



■nm 



-m 



22 21 

+ 3p + 4g)a + —m 2 p + —mpn — pn + —nm 4 + 5r — —m" 

5 5 5 5 



18 



2 2 
n m 



+ 2g H n m — 4m n H — m n + 3qm 

5 5 



3mr 



3 



'inq — 5gm 



2mp — -to 



12 o o 4 c 3 9 
—m i p + n A + —to - -n 2 
5 5 5 



= 0. 



(74) 
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Since p 32 is quadratic in a we solve it for a, and obtain, 

a = i(-13nm - 10n 2 + 4m 3 + 20q + I7m 2 n - 4m 4 + 15p - \7mp 

+ (-4Qqm A + 80qm 2 + 40m 3 p + 60np 2 - 15n 2 m 4 - WOmpn - 200nq 
- I5n 2 m 2 + 400g 2 + 60n 3 m 2 - 100n 2 q - 80mpn 2 + 200m 2 r + 260m 2 nq 
+ 225p 2 - 120m 3 r + 40m 5 p + 265mV - 40<?m 3 - 80m 4 p - 20gmn 
+ 360toVj + 30n 2 m 3 + 600p<2 - blOmp 2 - 120n 3 m - 680mpq 
+ mOmnr - 500nr + 80pn 2 - I70m 3 pn + 60n 3 ))i/(2m 2 - 5n). (75) 

So a is a number calculated directly from the coefficients of the quintic. Now 
the coefficient multiplying the linear term in d, p 3 i, is also linear in both r\ and 
£. Solving it for r\ and substituting this into the zeroth term in d, p 30 , gives a 
quadratic equation in £, i.e. let 

P30 = 6? 2 +£i£ + ?o=0 (76) 

and solve, giving, 

? = ^(-£i + & 2 -4fc£o)*)/6. (77) 

Where the equations for £ , £i an d £2, depend only on the coefficients of the 
quintic, and are given in the pre-print 27 . P2 is cubic in d, i.e. 

p 2 = d 3 d 3 + d 2 d 2 + did + d (78) 

solving using Cardano's rule on Maple, gives 

d = \{md 1 d 2 d 3 -lQ8d d 2 -8d 3 
6 

+ 12^3(4^3 - d\d\ - \8d x d 2 d 3 d Q 
+ 27 d%d 2 3 + 4d d 3 2 d 3 )i)i /da 

- f (3di<*3 - d%) I '(d 3 (36did 2 d 3 - 108d o d 3 - 8^ 
+ 12V3(4rf?d 3 - d\d\ - \8d x d 2 d 3 d {) 

+ 27 d 2 d 2 3 +4d d 3 d 3 )?)i) 

- \d 2 /d 3 . (79) 

Where do, d\, d 2l and d 3 , are also given in the pre-print 27 . Now (67) has 
been transformed to (72), where A and B depend only on the coefficients of 
the quintic(65). So that equation(72), with a linear transformation, becomes 
Bring's normal form, 

z 5 -z-s = (80) 



11 



where, 



(81) 



s = -B/{-A)i (82) 

We have used the fact that P4 is linear in a to get p 4 = 0. Then 6, c and d were 
considered a point in space, on the curve of intersection of a quadratic surface 
p 3 , and a cubic surface, p 2 . Giving p 3 = p 2 = 0, as required 24 ' 26 ' 28 . 



B. The Solution to Bring's Normal Form 

Bring's normal form is solvable. Any polynomial that can be transformed to 
the form, 

z n -az m -b = (83) 

can be solved with the hypergeometric equation 29 . The solution is given by 
considering z — z(s) and differentiating Bring's equation (80) w.r.t. s four 
times, then after making the substitution, 




(84) 



one obtains the fourth order Fuchian generalized hypergeometric differential 
equation, 

3/1 ,d 4 z 2 ,9 r7 ,d 3 z ,51 411 d 2 z ,3 183 .dz 231 
1 (1 - *> dF + t { 2 - ?t) ^ + t( T6 - 40 *> + ( 32 - lo *> * + 1^0000 Z = 



whose solution is 30 33 , 

z(s) = 4 Fs([ 20 , 35], [4, 4, 3]. 44O (86) 

C. Reverse Tschirnhausian Transformation 

When undoing the Tschirnhausian transformation one must be careful to gen- 
erate and keep all possible roots. All roots are calculated below. Incorrect 
roots are eliminated by direct substitution. Firstly, calculate y, with (81). We 
now undo the Tschirnhausian transformation by substituting d, c, 6, a and y 
into the quartic substitution, (64). The resulting equation is solved using the 
quadratic method 3 . The quadratic method introduces all the possible roots for 
the quartic. Let the general quartic, 

x 4 + a 3 x 3 + a 2 x 2 + aix + a = (87) 
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have coefficients, 

a = a + y 
a\ = b 
a 2 = c 
a 3 = d. 

(88) 

There are twelve possible roots to the quartic, let them be contained in the 
array, x mn with m € (1,2,3) and n e (1,2,3,4). Only four of these roots 
satisfy the quartic for a given set of coefficients ao, Oi, a 2 and a 3 . Using nested 
loops to enter the candidate solutions into both the general quintic (65) and 
the quartic(87), let the root that satisfies them both be r\. This root is then 
factored out of the quintic(65), giving the following equations, 

a = q + r\p + r\n + mr\ + r\ 
di = p + r\n + r\m + rf 
a,2 = n + mri + r\ 
a 3 = m + ri 

(89) 

These equations are then substituted into the general quartic, which gives an- 
other twelve roots, x mn . This time four of these roots satisfy both the quartic 
and quintic, let these roots be r 2 , r 3 , r± and r 5 . They were found the same way 
as ri, using nested loops. The final array has five non-zero elements, the five 
complex roots of the general quintic equation(65). 



D. Conclusion 

The five formulae for the roots of the general quintic equation are obtained, 
without root ambiguity, in agreement with Gauss' 'Fundamental Theorem of 
Algebra' when the polynomial is of degree five. 



VI. Determination of the Genus Two Hyperellip- 
tic Periods for Arbitrary Parameters 

Now from the hyperelliptic function(49), the periods can be calculated exactly 
for arbitrary branch places calculated in the previous section, this gives, 

W M = u k Cl,ai = phefk(ci,ai) + shef k (a,ai) (90) 

where 

(Ak,o + A k>1 ai) 



phefkicuai) = - 



iy/ (c - ai)(a 2 - ai)(c 2 - 01) 



( 5 \ 1 1 1 1 1 n C\ — ai Ci — cii C\ — cii 



x F D ^> (1, 1 -; 2; I^i, 1, 1) (9 1) 

2 2 2 2 2 c — ai s 2 -ai c 2 -«i 
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, , / \ A ktl (a - ai) 
8kef k (ci,ai) = 

2iyJ(c - ai)(a 2 - ai)(c 2 - ai) 

2 2 2 2 2 c — ai a 2 — a x c 2 — a x 

which becomes, 

, n(A kfi + A fc iai) (3) 1 111 

^^"^ F D (»>d, i.i.i^jx^,*,) (93) 



4^/(0 - ai)(a 2 - ai)(c 2 - ai) '" 2 2 2 2 
with, 

ci - ai ci - ai ci - ai . . 

.Ti = ;x 2 = ;x 3 = (94) 

a 2 — ai c 2 — ai c — ai 

and similarly, 

0J k , 2 = -Uk c ' a2 = -phef k {c, a 2 ) - shef k (c, a 2 ) (95) 

(A k ,o + A kA a 2 )Tr Fd (.S)A I I I-i-^ X2 X3 ) 

2\/( a i - a 2 ){c 2 - a 2 )(ci - a 2 ) 2'2'2'2' 
7rA M (c-a 2 ) (3 ) ,3 1 1 1 



and, 



W fc,2 = 



4\/ (°i ~~ a 2)(c 2 — a 2 )(ci - a 2 ) 2 2 2 2' 

c - a 2 c - a 2 c - a 2 

an = ;:r 2 = ;x 3 = (97) 

ai — a 2 c 2 — a 2 Ci — a 2 

w fc, 2 = Mfc C2;a2 = phef k (c 2 , a 2 ) + shef k (c 2 , a 2 ) (98) 

7r(A M +Afc,iQ 2 ) F D (3) ( X - - --1-xi x 2 x 3 ) 

2 V( a i ~ °2)( c + «2)(ci - o 2 ) 2 2 2 2 

it A k \{c 2 — a 2 ) „ t n i.3 1 1 1 „ , 

+ — 41 7 =F D (3) (-,-,-,-;2;an,x 2 ,x 3 ) (99) 

4v(ai - «2)(c - o 2 )(ci - 02 ) 2 2 2 2 

c 2 - a 2 c 2 - a 2 c 2 - a 2 

xi = ;z 2 = ;£3 = (100) 

ai — a 2 c - a 2 Ci — a 2 

and finally, 

Wfc.i = Ufc a,Cl = phef k (a, ci) + shef k (a, c\) (101) 
Wfc,i =Ufe C2,ai + Wfc,2 (102) 

where, 

2^( a 2 - «i)(ci - oi)(c - oi) 2' 2' 2' 2' 

+ — , = = F D { 3) (-, -,-,-; 2: an, x 2 , ar 3 ) (103) 

4V(a 2 - oi)(ci - oi)(c- ai) 2 2 2 2 
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x x = ;x 2 = ;x 3 = (104 

c-ai - 

Riemann's moduli are given by, 
T M = 

^2,2 = ^ — ^ (105) 

and by Riemann's condition(26), 

Tl.2 = 



C2 — ai 

x 2 = 

ci - ai 


C2 

;^3 = — 

«2 


<J->2,2U 1,1 — k 


'l,2W 2,1 


A 




— L02,l< J ->' 1,2 "+ 


- Wi^W 2,2 


A 




^2,iW 1,1 — a 


•'l, 1^2,1 


A 




— d- , 2,2 w 'l,2 "+ 


- Wi j2 W 2,2 


A 



= 75,1 (106) 

VII. The Solution to the Hyperelliptic Integral 

The solution to the hyperelliptic 4> integral(7)is with x\ = r g /r, 

r g 0i2 2 ^ai - a 2 



Cs/a[ - a 2 12 2 - yj{ai - ci)(ai - c 2 )(ai - c)9 13 2 



(107) 



where #i 3 and #12 are given in (60) and (62). The hyperelliptic <j> integral(10)may 
be written as 

f Xl Ai ixdx 

— = / — ; =phef(A ltl ;xi,a ) + shef(A 1>1 ;x 1 ,a ) (108) 

Z Jao V 

phef(A 1<1 ;x 1 ,a ) = B(a )A 1>1 x 

(l,mi +m 2 +m 3 + m 4 + m 5 )(|, mi)(i, m 2 )(|, m 3 )(|, mi)(\, m 5 ) 
,4"^ (2, mi + m2 + m 3 + m 4 + ra5)mi!m2!m3!m4!m5! 

(m)=0 

x X2 mi x 3 m2 x 4 m3 x 5 mi z m5+1 (109) 

s/ie/^i,!;^!,^) = ~ ao)A\,\ x 

(2, mi + m 2 +m 3 + m 4 + m h ){\, mi)(|, m 2 )(|, m 3 )(^, m 4 )(|, m 5 ) 
,•4^ (3, mi + m 2 + m 3 + m 4 + m 5 )mi!m2!m 3 !m 4 !m 5 ! 

(m)=0 

x x 42 mi x 3 " i2 x 4 m3 x 5 m4 ^ m5+2 (110) 

and 

4>' f Xl ( A 2,o + A 2 , 1 x)dx 

— = / = phef{A 2fi ,A 2! i;xi,a ) + shef(A 2>1 ;xi,a ) 

Z Ja a V 

(111) 
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phef(A 2 ,o, A 2> i;xi,ao) = B x (a A 2 . 1 + A 2fi ) 

(l,mi +m 2 +m 3 + m 4 + m 5 ){\, m{)(\, m 2 ){\, m 3 )(^, to 4 )(|, m 5 ) 
,"4^ (2, mi + m 2 + m 3 + m 4 + m 5 )mi!m 2 !m 3 !m 4 !m 5 ! 

(m)=0 

x x 2 mi x 3 m2 X4 m3 X5 m4 z m5+1 (112) 

shef(A 2 X;X 1 ,a ) = ^BA 2 . 1 (r il - a ) x 

(2, mi +m 2 +m 3 + m 4 + m 5 )(±, mi)(i, m 2 )(|, m 3 )(|, m 4 )(±, m 5 ) 
^ (3, mi + m 2 + m 3 + m 4 + m5)mi!m 2 !m 3 !m 4 !m5! 

(m)=0 

X £ 2 mi X 3 m2 X 4 m3 X5 m4 Z m5+2 (113) 

Abel's coefficients j4 2i o and A 2i i depend on the cosmological constant and are 
determined by experimentation and observation. So with 

ai - a ai-a a\ - a a\ - a 
«2 = ; £3 = ; £4 = ; 2^5 = (114) 

a 2 - a c - a Ci — a c 2 - a 

where the branch places are related to the apsidal points by, 

»i = — ;a 2 = -^;c = -^-;ci = —\c 2 = S- (115) 
n r 2 r 3 r 4 r 5 



and the variable place is given by 

xi - a 



(116) 



ai - a 
with 

|sb 2 | < 1; \x 3 \ < 1; M < 1; M < 1; \z\ < 1 (117) 
so that the variables in the argument of the theta functions are, 

— LOl 2 6' + L0 2 2 (j) 

Vl 55 ' A ' 

v 2 ee "i.^ (us) 

where the periods u>ij are given in section VI. With this choice the equations of 
Forsyth and Baker arc identical to these, in the case of g = 2. The roots of the 
quintic ai, a 2 , c\, c 2 and c are given by setting the coefficients in(41) equal to 

4Am 2 c 2 r„ 4 

A ° = 4r = 3L 2 
Ai = 4q = 



4 A 2 4m 2 c 2 r 2 4r 2 .E 2 
4m 2 c 2 r 2 



A 3 = 4n = 



L 2 

A 4 = 4m = -4 (119) 
then substituting into the solution to the quintic in section V. 
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A. Discussion 



The geodesies of a particle in SdS space are determined by three coordinates, 
the radial coordinate, r, and two angular coordinates (p and <j>' . <f> is the actual 
physical coordinate, while </>' is an angular coordinate that comes from Abel's 
theorem, r is a function of both <fi and <f>', by the quadrupoly periodic thcta 
functions of two variables. Now <f> and <f>' depend on r, i.e. <j> = 4>(r) and 
(jjf — (j)'(r) and they can be expressed as a linear combination of phef and shef. 
Thus (j) = <p(r) = —2phef(r) — 2shef(r) is the closed form solution to the 
calculation of the geodesies in SdS space. 

VIII. The Reduction of the Hyperelliptic Func- 
tions to Elliptic Functions when A = 

When the cosmological constant is set to zero and the particles are subject to 
gravitation alone, two roots of the quintic(42) go to zero. Let these two roots be 
c 2 and a 2 . Then iv'1.2 = oj' 2 , 2 = 0. Then the following equations are obtained 
for the Riemann moduli, 

Tl,l - (W2,2w'l,l - wi, 2 w' 2 ,i)/A ^ (120) 

T2,2 = (-^2,l w 'l,2 + Wl,iw'2,2)/A = (121) 
n,2 = (^2,l w 'l,l - Wl,l^'2,l)/ A = (-W2,2w'l,2 + W1,2W' 2 ,2)/A = (122) 

Let Abel's coefficients be functions of the cosmological constant, i.e. 

A 1A = Ai,i(A), A 2fi = A 2 , (A), A 2il = A 2A {A) (123) 
Taking the limit as A -» with (108), (111), (118), one has 

-U 12 <t>' + ^2.20 „ n „.s 

vi = lim : : — = 124 



gives, 



and 



or, 



A 2 , (0) = (125) 
A 2 ,i(0) = ^£ (126) 

Wl,2 

= Um ^'- ^ = Um ±_ (12?) 
A^O A A^O LJi.2 

_2 

v 2 = lim (phef (A hl (A);x,y) + shef (A ltl (A);x,y)) (128) 

A^O Uli 2 
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Letting liniA^o ^i,i(A) = Ai^(0) = Wi^ and using the identity, 

ZF D ] (! . \ \, \, \\ : 2 ; X 2 , X, X, X 3 , Z) 
1 n (5), 1 1 1 1 1 

--.tz (2, -,-,-,-,-; 3; x 2 ,x,x,x 3 ,z) 

= zF$\\,\,\±2-,x 2 ,x z ,z) (129) 
to obtain the tripple hypergeometric function, 

V2 2i v /( ai -y)( Cl - 2/ )(c-y) D 1 '2'2'2' ' Cl - y ' c - y ' oi - y ' 

(130) 

that reduces to the incomplete elliptical integral of the first kind 35 36 , when 
y — * ai. So that becomes, 



v 2 = 4> = 



With these definitions, 



F { sin-\[*EK, X [^K) (131) 
V ci — a\ V c — ai 



where sn(u) is Jacobi's elliptic function 37 . Now with the choice 



c = —,ci = —,ai = —,k= * , 133 

r\ r 2 r 3 V TiT\ — r 3 r 2 

which gives the Schwarzschild black hole solution without the cosmological con- 
stant, as a genus one elliptic function, which is the same as Ref. 7, it is 

(---) = (--->n 2 M (134) 
r n r 3 n 

Lastly, the moduli become 1 and 
w'2,1 



T i,i 



_ limA^ (p/ie/ 2 (ci,ai) + she/2 (ci, ai)) 



liniA^o^ 2 " 1 +^ 2 , 2 ) 

Using the identities, 



(135) 



Fh{\W\-^x,x,z)- -xFl{~,~;2;x,x,z)= ^(H;^) (136) 
*F<g\\, \, \, \; 1; 1, x, y) - tt^ 3 (| ± ± ±; 2; 1, x, y) = -| F(sm -i^ ^lST) 
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where 2^1(5, 5; 1; z) is the Gauss function for the complete elliptical integral of 
the first kind and F is the elliptic function for the incomplete elliptical integral 
of the first kind. Then the periods of section VI, (90) to (104)reduce to elliptic 
functions, i.e., 

' a 1 c^ 2 )/ 1 1 1 3 1 c i _a i\ 

y/c — a\ 2 2 2 2 c — a\ 

1 tt ,11 
= A fe ,i— ===- 2 Fi(-,-;l;z) 
y c - fli z 2 2 



= A k ,^J=F(* x ^^), (138) 



similarly, 

C2,ai 0,ai 



Ak,i—i=^=F{sin~ 1 ^j—^^—, \ — — — ) (139) 



Wfc , 2 = -A fc . 1 ^=^( S m- 1 (./^— ),,/^i) (140) 
Vc — ci V c — ai V c — c i 

=w°' Ql +^,2 (141) 

and lastly, 

w'fc,2 - (142) 
It has been shown that the solution for the geodesies in SdS space in terms of 
multiple hypergeometric and hyperelliptic modular functions reduce to elliptic 
functions in the case of pure Schwarzschild space. 



IX. Conclusion 

The geodesies of particles in the combined Schwarzschild black hole and a non- 
zero cosmological constant, are given by ratios of quadrupoly periodic thcta 
functions of two variables. The inverse functions are obtained as multiple hy- 
pergeometric functions of five variables, thus solving the Inversion Problem for 
hyperelliptic integrals of the first kind, i.e. genus two. It is a closed form so- 
lution since the branch places are found by solving the quintic equation for it's 
five roots. The five formulae for the five roots of the quintic were obtained, in 
agreement with Gauss' Fundamental Theorem of Algebra when the polynomial 
is of degree five. The closed form solution for the geodesies in SdS space was 
found to be <fi = 4>{r), which is a linear combination of phef(r) and shef(r) 
which are convergent twelve slotted five variable hypergeometric functions. 
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